The band structure of hexagonal diborides ZrB 2 , VB 2 , NbB 2 and TaB 2 
in comparison with the superconducting MgB 2 
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The band structure and the Fermi surface of hexagonal diborides ZrB2, VB2, NbB2, TaB2 have 
been studied by the self-consistent full-potential LMTO method and compared with those for the 
isostructural superconductor MgB2 . Factors responsible for the superconducting properties of AIB2- 
like diborides are analyzed, and the results obtained are compared with previous calculations and 
available experimental data. 
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A recent discovery [1] of the critical transition (T c sa 
40 K) in magnesium diboride (MgB2) and development 
of some promising superconducting materials on its ba- 
sis (ceramics, films, extended wires, see review [2]) gave 
an impetus to active search for novel superconductors 
(SC) among related compounds with the structure and 
chemistry similar to MgB2. It is reasonable that hexago- 
nal (AlB2-like) metal diborides isostructural with MgB2 
were considered as first SC candidates. Detailed inves- 
tigations of the band structure and coupling mechanism 
in MgB 2 [2-7] showed that diborides of group I, II met- 
als of the Periodic System, for instance metastable CaB2 
[6], LiB 2 , ZnB 2 [7], hold the greatest promise as SC can- 
didates. The authors of [8] predicted the possibility of 
critical transition with T c > 50 K in AgB2 and AuB2. 
It is much less probable to find new SC (with T c > 1 
K) among AlB2-like diborides of d metals (MB2), see [2- 
7]. The first report ([9], 1970) on the superconductivity 
in NbB2 (T c ks 3.9 K) was not confirmed by systematic 
studies [10] of SC properties of a series of diborides MgB2 
(M = Ti, Zr, Hf, V, Nb, Ta, Cr), according to which their 
T c are smaller than 0.7 K. 

Therefore recent reports [11-13] on rather high T c for 
ZrB 2 (5.5 K [11]), TaB 2 (9.5 K [12]) and NbB 2 (5.2 K 
[13]) were quite unexpected. It is remarkable that in the 
investigations of identical series of diborides (TiB 2 , ZrB 2 , 
HfB 2 , VB 2 , NbB 2 , TaB 2 [12] and ZrB 2 , NbB 2 , TaB 2 [11]) 
each group revealed "its" superconductor, namely ZrB 2 
[11] and TaB2 [12], while all other MB2 phases were as- 
signed to non-superconductors. 

The results published in [12] prompted the authors of 
[14] to examine in detail the temperature dependencies of 
magnetic susceptibility and electrical resistance in TaB2. 
It was established that the SC transition for TaB2 is not 
observed to T w 1.5 K. The SC properties of MgB 2 and 
TaB2 discussed [14, 15] on the basis of the band structure 
calculations were found to be essentially different due 
to strong hybridization effects of the Ta5d-B2p states. 
A weak (as compared with MgB2) electronic interaction 
with the E2<7 mode of the phonon spectrum was noted 
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in [14]. An abrupt lowering of T c for TaB2 (and the 
absence of SC for VB2) is explained in [15] by a consid- 
erable decrease in the contributions from B2p states to 
the density of states at the Fermi level (N(Ejr)): MgB 2 
(0.494) > TaB 2 (0.114) > VB 2 (0.043 states/eV). Ana- 
lyzing fine features of soft X-ray emission and absorption 
BK-spectra of MgB2, NbB2 and TaB2, the authors of 
[16] pointed out fundamental differences in the structure 
of their pre-Fermi edges with dominating contributions 
from B2p CT (MgB 2 ) or B2p T states (NbB 2 , TaB 2 ). Any 
operations reproducing results [11], to us are unknown. 
As is shown in [2-7] , the superconductivity in MgB2 and 
related borides is fairly well described in the context of 
the electron-phonon interaction theory. Therefore the pe- 
culiarities of the electronic spectrum, primarily the com- 
position and structure of the near-Fermi bands, consti- 
tute the major factor responsible for the formation of this 
effect. In this work we present the results of detailed band 
structure analysis of Zr, V, Nb and Ta diborides in com- 
parison with the superconducting MgB2. As is known, 
these diborides are isostructural (A1B 2 type, space group 
P6/mmm), their crystal lattices are made up of alternat- 
ing hexagonal metal (M) monolayers and graphite-like 
boron layers [17]. The unit cell contains three atoms (M, 
2B). The main differences are due to the metal sublatticc 
type, viz. the electronic configurations of metal atoms 
(Mg - 3s 2 3p°; Zr - 5s 2 4d 2 ; V, Nb, Ta - (n+l)s 2 nd 3 , where 
n = 3, 4, 5 respectively), which determine the electron 
concentration (EC) growth (MgB 2 (8) < ZrB 2 (10) < 
VB2, NbB2, TaB2 (11 e/cell)) and the variation in inter- 
atomic bonds, see [2, 18, 19]. 

The band structure of MgB 2 was calculated in the frame- 
work of the LDA theory by the self-consistent full- 
potential linear muffin-tin orbital (FP LMTO) method 
with allowance for relativistic effects and spin-orbital in- 
teractions [20, 21] with the exchange-correlation poten- 
tial in accordance with [22]. Equilibrium parameters of 
MB2 cells (Table 1) were obtained based on the condition 
of the system's total energy minimum. 
MgB2. As follows from Figs. 1 and 2, the peculiarities 
of the band structure of the superconducting MgB2 are 
due to the B2p states, that form four a(2p x ^ y ) and two 
7r(p 2 ) energy bands. The B2p z states are perpendicular 



2 



TABLE I: Diboride lattice parameters for Mg, Zr, V, Nb, Ta 
by our FP LMTO self-consistent calculations and experiment 
results 



Diboride 


; a (A) 


c (A) 


c/a 


a* (A) 


c* (A) 


c/ a* 


MgB 2 


3.04869 


3.46637 


1.1370 


3.083 


3.521 


1.142 [23] 


ZrB 2 


3.16932 


3.53126 


1.1142 


3.170 


3.532 


1.114 [11] 










3.165 


3.547 


1.120 [23] 


VB 2 


3.00678 


3.04768 


1.0136 


2.997 


3.056 


1.0196 [15] 










2.998 


3.057 


1.020 [23] 


NbB 2 


3.18141 


3.35693 


1.0550 


3.116 


3.264 


1.060 [23] 


TaB 2 


3.16421 


3.32337 


1.0503 


3.082 


3.243 


1.0522 [12] 










3.098 


3.224 


1.0407 [14] 










3.083 


3.244 


1.0522 [15] 










3.097 


3.225 


1.041 [23] 



- [11,12,14,15,23] 



to boron atom layers and form weak interlayer ir bonds. 
The B2-p x y bands are of two-dimensional-like (2D) type 
and form flat areas in the T-A direction of the Brillouin 
Zone (BZ). A small dispersion of a bands is also indica- 
tive of an insignificant interaction between Mg-B layers. 
Two B2pj;. y bands intersect and make an appreciable 
contribution to the density of states (DOS) at the Fermi 
level being responsible for the metal-like properties of 
MgB 2 ,see Table 2. One of the most important features 
of MgB 2 is the presence of hole B2p :Cjy states: in the T-A 
direction they are above E p and form cylindric hole- type 
elements of the Fermi surface (FS), Fig. 1. 

Thus, the distinguishing characteristics of the band 
spectrum of MgB2, that are crucial for its superconduct- 
ing properties, as well as inter- and interlayer interaction 
effects (see also [2-7]), include: (1) location of cr(p x , y ) 
bands relative to E F (the presence of hole states); (2) the 
value of their dispersity in the T-A direction (AE CT (r-A), 
determined by the degree of interaction between metal- 
boron layers); (3) the value and orbital composition of 
N(E^) (dominating contribution of s states of boron 
atoms from graphite- like layers). Let us consider in this 
context the band structure of Zr, V, Nb and Ta diborides. 
First note that the most obvious consequence of the vari- 
ation in the composition of the metal sublattice in the 
series of diborides is an increase in the EC with succes- 
sive occupation of energy bands. Then the Fermi level 
for ZrB2 is located in the pseudogap between completely 
occupied bonding and free antibonding states. This de- 
termines the maximum stability of ZrB2 (and also of iso- 
electronic and isostructural TiB 2 and HfB 2 ) in the series 
of AlB2-like phases and their extreme thermomechanical 
characteristics [23]. These inferences were confirmed by 
recent FP LMTO calculations of the cohesion energy of 
some MB 2 phases (M = 3d, 4d and 5d metals) [18, 19]. 
MgB 2 and ZrB 2 . As is seen from Figs. 1, 2 and Table 
2, the structure of the near- Fermi spectra edges of ZrB 2 
and the SC MgB 2 differ radically. For ZrB 2 , (1) cr(p x , y ) 
bands of boron are located in the region are located below 
Fjp (-1-1 eV in point A BZ) and the corresponding hole 



TABLE II: Dispersion of a(P x , y ) bands in the direction T-A 
(AE CT (P-A), eV) and orbital contributions to the density of 
states at the Fermi level (states/eV, cell) in A1B 2 phases. 





AE CT (r-A), 








DOS 








Phase 


eV 


Total 


Ms 


Mp 


Md 


Mf 


Bs 


Bp 


MgB 2 


0.72 


0.719 


0.040 


0.083 


0.138 




0.007 


0.448 


ZrB 2 


1.73 


0.163 


0.001 


0.002 


0.130 




0.000 


0.030 


VB 2 


2.69 


1.379 0.024 0.013 


1.255 




0.002 


0.085 


NbB 2 


2.49 


1.074 


0.037 


0.017 


0.818 




0.012 


0.190 


TaB 2 


2.61 


0.910 


0.003 


0.016 


0.664 


0.038 


0.011 


0.178 



states are absent; (2) a considerable dispersion of these 
bands appears in the T-A direction ((AE a (r-A) = 1.73 
eV) , cr bands become no longer of the 2D type as a result 
of formation of strong covalent d-p bands between metal- 
boron layers, in which partially occupied 7r(p z ) bands 
take part; (3) the value of N(E^) decreases drastically 
as compared with MgB2 (from 0.719 to 0.163 states/eV), 
the maximum contribution (~80%) to N(E^) being made 
by Zr4d states (whereas the contributions from boron 
states are much smaller - ^18%). The change in the type 
(2D — > 3D) of the near-Fermi states can be easily traced 
by comparing the structure of the FS of MgB2 and ZrB2, 
Fig. 1. It is seen that the FS of ZrB 2 consists of three 
types of figures defined by mixed Zr4d,5p-Bp states: (a) 
a 3D rotation figure around a straight line along the T-A 
direction with the hole-type conductivity; (b) a 3D figure 
near the centre of the M-K segment with the electronic- 
type conductivity; and (c) tiny 3D-type sections with the 
electronic conductivity. 

VB 2 , NbB 2 and TaB 2 . Energy bands, FS and DOS of 
these isoelectronic and isostructural diborides are demon- 
strated in Figs. 3, 4, and some electronic structure 
parameters are listed in Tables 2, 3. The mentioned 
above differences in ZrB 2 and the SC MgB 2 (filling of 
<j(\>x,y) bands, decreased contributions of B2p states to 
N(Ep), variation in the type (2D — > 3D) of the near- 
Fermi states) are typical also of VB 2 , NbB 2 , TaB 2 . Be- 
sides, they have the following features in common (as 
compared with ZrB2): (1) partial occupation of the an- 
tibonding d band responsible for the metallic-type con- 
ductivity; (2) considerable growth of N(Ep); and (3) in- 
creased filling of 7r(p z ) bands. A peculiar form of the 
surface Fermi transformation is observed for example for 
TaB2 (Fig. 3): the surface Fermi contains double ("in- 
ternal" and "external") non-intersecting electronic-type 
rotation spheroids around point A defined by 3D-B2p 
and Ta5d xz . yz states respectively. The bonding a(p XtV ) 
bands of boron lie in the region of -1.3, -2.5 and -2.6 
eV in A-point of BZ for VB 2 , NbB 2 and TaB 2 accord- 
ingly, below Ep and have, as in the case of ZrB2, a sig- 
nificant energy dispersion AE CT (r-A), which has a maxi- 
mum value for VB 2 , Table 2. In the series of isoelectronic 
VB 2 -> NbB 2 -> TaB 2 , N(E F ) decreases systematically 
and has a maximum value for VB2 due to the contribu- 
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TABLE III: Density of states at the Fermi level (states/eV) 
for Mg, Ta, V diborides according to our FP LMTO calcula- 
tions in comparison with the results [14, 15]. 







N F 




Phase 


Our data 


FP[15] 


FPLO[14] 


MgB 2 


0.719 


0.691 


0.71 


B2p* 


0.448 


0.494 




TaB 2 


0.91 


0.966 


0.91 


Ta5d- 


0.664 


0.647 




B2p- 


0.178 


0.114 




VB 2 


1.379 


1.359 




V3d- 


1.255 


1.235 




B2p- 


0.085 


0.043 





* - total DOS and partial contributions of Md and B2p 
states 



tion from the near- Fermi quasi-plane Y2>A XZtVZ band in 
the direction T - A. By contrast, the contribution of B2p 
states (antibonding s and p bands) to N(E^ ) in this series 
changes non-monotonously and reaches its peak (0.190) 
for NbB2, which is much smaller than that for MgB2 
(0.448 states/eV). A great concentration of B2p states in 
the vicinity of for NbB 2 (as compared with TaB 2 ) is 
also evident from spectroscopy experiments [16]. 
Thus, the performed analysis of the band structure and 



the surface Fermi of isostructural d-metal (Zr, V, Nb, Ta) 
diborides allows us to formulate their fundamental differ- 
ences from those of the superconducting MgB2: (1) oc- 
cupation of bonding p^ bands and the absence of hole s 
states; (2) increased covalent interactions between boron 
and metal layers (due to the hybridization of B2p-Md 
states) and the loss of the energy bands two-dimensional 
character; (3) changes in the value and orbital compo- 
sition of N(Ej?), where the dominant role is played by 
valence d states of metals. The latter circumstance is 
typical of low-temperature superconductors, for example 
metal-like compounds of these d elements with carbon, 
nitrogen, silicon (NbN, VsSi etc.), the T c values of which 
correlate with the values of N(Ej?) [24]. In this case the 
results obtained suggest that V, Nb, Ta diborides are 
more likely to have low-temperature superconductivity, 
and among them the maximum T c value can be antici- 
pated for VB2. On the contrary, if we assume that the 
major electronic factor for the formation of SC proper- 
ties of MB2 is the near-Fermi density of B2p states (by 
analogy with MgB 2 [2-7]), NbB 2 should possess the high- 
est critical transition temperature. It is noteworthy that 
according to the coupling model proposed in [25, 26], 
not only boron s, but also p bands occupation should 
be considered in this case. Anyhow the superconducting 
transition for ZrB2 is the least probable, and the results 
[11] need to be revised. 
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